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Abstract 

Two-dimensional turbulence governed by the so-called a turbulence equations, which 
include the surface quasi-geostrophic equation (a = 1), the Navier-Stokes system 
(a = 2), and the governing equation for a shallow flow on a rotating domain driven 
by a uniform internal heating (a = 3), is studied here in both the unbounded and 
doubly periodic domains. This family of equations conserves two inviscid invariants 
(energy and enstrophy in the Navier-Stokes case), the dynamics of which are be- 
lieved to undergo a dual cascade. It is shown that an inverse cascade can exist in 
the absence of a direct cascade and that the latter is possible only when the inverse 
transfer rate of the inverse-cascading quantity approaches its own injection rate. 
Constraints on the spectral exponents in the wavenumber ranges lower and higher 
than the injection range are derived. For Navier-Stokes turbulence with moderate 
Reynolds numbers, the realization of an inverse energy cascade in the complete 
absence of a direct enstrophy cascade is confirmed by numerical simulations. 
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1 Introduction 



Incompressible fluid turbulence in two dimensions has always been a subject 
of interest because of its relative simplicity and broad applications in geo- 
physical problems. Although the theory of two-dimensional (2D) turbulence 
is not strictly applicable to a geophysical fluid, the third (vertical) dimension 
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in many geophysical systems can nevertheless be "removed", and the result- 
ing 2D equations are thought to be quantitatively adequate for describing the 
large-scale horizontal motion. 

Some systems in this category, which receive considerable attention, are the 
surface quasi-geostrophic (SQG) equation, describing the motion of a rotating 
stratified fiuid [3,15]; the "rotating shallow flow" (RSF) equation, govern- 
ing the motion of a thin fluid on a rotating domain [35]; and the Charney- 
Hasegawa-Mima equation, governing the potential vorticity of an equivalent- 
barotropic fluid [6,13,14]. Except possibly for the last equation, these systems, 
together with the more familiar 2D Navier-Stokes (NS) equations, can be 
studied within a unifled framework by a family of models (known as the a 
turbulence equations) indexed by a parameter a e (0, oo) [23]. Like the 2D 
NS system, the nonlinear transfer of a turbulence conserves two quadratic in- 
variants: *a = /o°°*a(^)dfc and *2a = /o°° *2a(^) d/c, where = k^^{k), 
^'(/c) is the streamfunction power spectrum, ^ is a real number, and k is the 
wavenumber. Note that ^2(^) is the kinetic energy spectrum, and ^2 is the 
kinetic energy. 

The statistical dynamics of unbounded 2D NS turbulence are beheved to 
achieve a quasi-steady state which simultaneously exhibits a direct enstrophy 
cascade to high wavenumbers and an inverse energy cascade to low wavenum- 
bers. According to the classical theory formulated by Kraichnan [16,17], Leith 
[18], and Batchelor [2] (hereafter referred to as KLB), if the system is excited 
around a (forcing) wavenumber s, the dual cascade is supposed to carry virtu- 
ally all of the energy input to ever-lower wavenumbers, down to wavenumber 
A; = 0, and the enstrophy input to a high wavenumber around which 

the enstrophy is dissipated. Analyses based on Kolmogorov's phenomenology 
predict that the energy spectrum scales respectively as k'^^^ and k~^ in the 
energy- and enstrophy-cascading ranges (also known as inertial ranges). The 
KLB theory, when applied to a turbulence, implies the possibility that 
cascades to low wavenumbers and that cascades to high wavenumbers. 
These (possible cascading) invariants correspond to different physical quanti- 
ties for different ffows. For SQG turbulence, the direct-cascading candidate is 
the energy. For RSF turbulence, the direct-cascading candidate is the palin- 
strophy. When a = 4, the enstrophy becomes the inverse-cascading candidate. 
NS turbulence features the energy and enstrophy as the possible inverse- and 
direct-cascading quantities, respectively. 

Despite numerical evidence confirming the realization of an inverse energy 
cascade and of the Kolmogorov-Kraichnan k'^l"^ spectrum [4,5,12,19,32], the 
dual cascade remains conjectural. No convincing evidence exists of a direct 
enstrophy cascade. Rather, there exist some negative results with respect to 
this problem. In particular, for a bounded system in equilibrium, a direct en- 
strophy cascade is excluded [33]. For an unbounded system, a direct enstrophy 
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cascade is prohibited for all weak inverse energy cascades [32] . By "weak" , it is 
understood that the inverse cascades do not carry virtually all of the energy in- 
put to the large scales. This behaviour is presumably the dynamical behaviour 
for finite- Reynolds number turbulence. The questions then arise whether the 
inverse energy cascade can become "strong" in the limit of infinite Reynolds 
number and whether a direct enstrophy cascade can then be realizable. 



This study generalizes various previous results for the 2D NS system [30,32,33] 
to a turbulence and re-examines the dual-cascade hypothesis under the general 
framework of the latter case. Also, some new results on related issues, such as 
the non-conservative transfer of quadratic quantities other than the invariants, 
are reported. First, the derivations of some a priori estimates for the dissipa- 
tion terms of the quadratic invariants are presented. Second, the result of Tran 
and Shepherd [33] is generalized to the present problem, showing that under 
the infiuence of a viscous operator, the invariants satisfy < s^^a- Third, 
the conservative transfer of the quadratic invariants and the non-conservative 
transfer of other quadratic quantities are examined. The former consideration 
is essentially a reformulation of the dual-cascade hypothesis. It is shown that 
even in the limit of infinite Reynolds number, the inverse cascade of '^^ could 
become strong, carrying virtually all of the injection of to ever-larger scales, 
in the complete absence of a direct cascade of ^la- This study does not attempt 
an alternative to the classical dual-cascade theory, but is concerned with the 
plausibility of the dual-cascade hypothesis. A constraint on the spectral slopes 
for fluids in equilibrium is obtained by extending the derivations of Tran and 
Bowman [30,32] to the present case. Finally, the results from numerical simu- 
lations of NS turbulence are presented, confirming the realization of an inverse 
energy cascade in the absence of a direct enstrophy cascade. These numerical 
results suggest that for high Reynolds numbers, the turbulence would become 
non-direct-cascading by evolving through a transient direct-cascading stage. 



The remainder of this paper is organized as follows. Section 2 describes the 
equation of motion governing a turbulence and its conservation laws. Section 3 
presents the derivations of the upper bounds on the asymptotic averages of the 
dissipation of quadratic invariants and the dynamical constraint '^2a < s"\E'q,. 
Section 4 reports an analysis concerning the plausibility of a direct cascade in 
unbounded fluids. Section 5 features a proof of the creation and depletion of 
quadratic quantities other than the invariants. Section 6 presents a constraint 
on the spectral slopes of the kinetic energy spectrum for fluids in equilibrium. 
Section 7 reports the results from numerical simulations of 2D NS turbulence, 
conflrming the realization of a weak inverse energy cascade in the complete 
absence of a direct enstrophy cascade. The paper ends with some concluding 
remarks in the final section. 



3 



2 Governing equation 



The term "a turbulence" refers to a general class of 2D incompressible fluid 
turbulence, in which the conserved scalar (— A)°/^'0 is advected in a velocity 
field with stream function Here, a is a positive number, and A is the 2D 
Laplacian. The equation governing the fluid motion is 



9*(-A)"/^V' + Ai^, (-A)"/'V') = -v,{-/\Y{-/\r/^i^ + /, (1) 

where J(</?, (p) = d^fdycf) — dxfpdyip, z/^ > is the generalized viscosity coeffi- 
cient, > is the degree of viscosity, and / is the forcing. The stream function 
ip{x,y,t) is related to the fluid velocity u hy u — {—dyip,dxip). This family 
includes the three physically realizable fluid systems, which are the focus of 
this work. The case a = 1 corresponds to the SQG equation from the geophys- 
ical context, describing the motion of a stratifled fluid with small Rossby and 
Ekman numbers [3,15,23]. The familiar 2D NS system corresponds to a = 2. 
The case a = 3 corresponds to the RSF system, governing a shallow flow on 
a rotating domain with uniform internal heating [35]. 

The nonlinear term in (1) obeys the conservation laws 



{^J{^, i-Ar^'^)) = ((-A)"/VJ(V^, (-A)"/^)) = 0, (2) 

where (•) denotes the spatial average. As a consequence, the two quadratic 
quantities ^„ = KK-A)"/^^^) and *2a = i(l(-A)°/^V'P> are conserved by 
nonlinear transfer. 

In this study, / is assumed to be spectrally confined to a finite band of the 
wavenumber K = [ki, /c2], where ki > 0. More restrictive conditions on / and 
specific values of fi for the particular cases to be considered will be stated in 
due course. In particular, fi = 1 represents molecular viscosity as in the NS 
case. The SQG system has n = 1/2 as the physical dissipation mechanism; 
other values of /j, have also been considered for numerical purposes in the 
literature [21,25,26,28]. It is not well understood how the dynamics would 
depend on the degree of viscosity fi. Nevertheless, there appears to be a belief 
that z/^(— A)^ (for large positive values of fi) globally acts on the small scales, 
hence the name "small-scale dissipation". This belief is associated with the 
concepts of direct cascade and of dissipation range and is incorrect if a direct 
cascade is not realizable. For example, for 2D NS turbulence in equilibrium, 
with A)^ (/i > 0) the only form of dissipation, the dissipation of energy 
(enstrophy) occurs mainly in the large scales (around the forcing scale) [33]. 
This example means that there exists no dissipation range well separated from 
the forcing scale and that the viscous operator i/^(— A)'' would have to be 
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correctly termed "large-scale dissipation". In this case, the degree of viscosity 
II plays an important role in the spectral distribution of energy [30,33]. 

The problem of nonlinear transfer and spectral distribution of quadratic quan- 
tities such as energy and enstrophy is best handled in wavenumber space. 
Hence, it is convenient to formulate the present problem accordingly. The 
power spectrum "^^{k), which represents the power density of (—Ay^^ip (for 
real 9) associated with the wavenumber k, is defined by 

^e{k)^l I \{-Af'm\'de^lk' I m\'d£, (3) 

\e.\=k \£\=k 

where ip{£) is the Fourier transform of ip, and the integral is over all wavevec- 
tors £ having magnitude k. In this paper, most results apply to both un- 
bounded and periodic (bounded) cases; the integral over continuous wavenum- 
bers is then understood as the corresponding sum over discrete wavenumbers 
for the bounded case. The evolution of ^a{k) is governed by 

^^^a{k) = T^{k) - 2u^e'^<i!^{k) + F{k). (4) 

Here, T„(A;) and F{k) are, respectively, the spectral transfer and input func- 
tions. (See [11,16] for the famihar NS case.) By virtue of the conservation laws 
(2), the transfer function Ta{k) satisfies 

oo oo 

J (k) dk^J k'^T^ (k) dk = 0. (5) 





3 Asymptotic analysis 

In this section, several a priori estimates of some relevant dynamical quan- 
tities are derived, and their implications for power-law scaling of the energy 
spectrum are discussed. The results obtained include upper bounds on the 
asymptotic averages of the dissipation of a-^d a-^d an upper bound on 
^^20 in terms of and the forcing wavenumber. 

On multiplying (1) by ip and (— A)"/^-?/; and taking the spatial averages of 
the resulting equations, noting from the conservation laws that the nonlinear 
terms identically vanish, one obtains evolution equations for a-nd "if 2a, 
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^^« = -2l/^^r2^+„ + (^/), (6) 
= -2v,'^2,^2a + ((- A)'^/^^/) ■ (7) 

Using the Cauchy-Schwarz and Young inequalities and because the forcing 
spectral support is confined to the wavenumber interval [/ci,/c2], one obtains 
the upper bounds on the injection terms in (6) and (7): 



(V/) < (|(-A)'^/2W4^|2^1/2^|(_^)-/./2-a/4^|2^1/2 

<^.^2,+.+^(i(-A)--/2-"/vr) 

S i^/i"±'2M+Q + 72^' 

^'^/i" 1 

((-A)"/2^/) < (|(-A)''/2W2^|2^1/2^|(_^)-/./2_y|2^1/2 
<^.*2,+2a + ^(|(-A)-'^/Vr) 

<^M*2,+2a+^^. (8) 

The substitution of (8) in (6) and (7) yields 



d (IfP) 
These inequalities readily admit the upper bounds on the asymptotic averages 

of *2M+a and *2M+2a: 



lim 



sup i / dr < ^-Jp, lim=up 1 /(l/R dx, (10) 

M ^ 

t * 

If 1 1 

hm sup - / ^2M+2a dr < ^ hm sup - 

^ 



(l/r)dr. (11) 



For a time-independent forcing /, the asymptotic average on the right-hand 
sides of (10) and (11) reduces to 

The upper bound (10) on the asymptotic average of \l/2/i+a is interesting, par- 
ticularly for the unbounded case. This bound prohibits the spectrum ^'2/i+a(^) 
[\E'a(/c)] from acquiring a spectral scaling not shallower than [k"^'^'^] in 
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the limit k ^ because otherwise, a divergence of would entail. This 

means that a constant inverse flux of \E'o, (if it is realizable) ought to pro- 
ceed toward A; = (in the limit t — > oo) via a spectrum ^„(/i;) shallower than 
^-2^-1 'pjjjg observation is intuitively obvious from the physical point of view: 
the proceeding of a sustainable inverse cascade is incompatible with a persis- 
tent increase in the dissipation of the inverse-cascading quantity. Note that 
the critical scaling ^^(/i;) oc /c~^^~^ (for k < ki) may be realizable only down 
to a finite wavenumber ko > and necessarily corresponds to an asymptotic 
absolute equilibrium of 

For the NS system {a ~ 2 and fj, = 1), ^2^l+a = ^4 (the enstrophy) is the 
dissipation agent of the energy "$2- The above analysis implies that the en- 
strophy (energy) spectrum in the inverse- cascading region is no steeper than 
k'^ {k~^). This constraint turns out to be "generous", as the Kolmogorov- 
Kraichnan k~^^^ energy spectrum is much shallower than k~^. For SQG tur- 
bulence {a — 1 and = 1/2), the energy ^2 is the dissipation agent of 
the inverse-cascading candidate ^1. Hence, the energy asymptotic average is 
bounded. This condition requires that ^'2(^) in the low- wavenumber region 
be no steeper than k~^. Tran and Bowman [31] have previously noted this 
requirement. 

Remark 1. Like the constraint in the NS case, the constraint on the spectral 
steepness of the low-wavenumber region for SQG turbulence is the thresh- 
old at which an inverse cascade is not sustainable. To allow for a steady 
positive growth rate of ^1, an energy spectrum shallower than k"^ in the 
low-wavenumber range is required. 

A result derived by Tran and Shepherd [33] for NS turbulence confined to a 
doubly periodic domain is now generalized to the present problem. For sim- 
plicity, it is assumed that the forcing satisfies s°'{il)f) = ((— A)"/^'?/'/), where 
s is some (constant) wavenumber. (The result to follow can be readily gener- 
alized to a less restrictive forcing class [33].) This condition can be realized in 
several ways. The simplest example [33] is the case of a monoscale forcing at 
a wavenumber s, i.e. a forcing / for which (— A)^/ = s^^f. Another example 
[27,30,32,33] is 



ioT k & K — [ki, /C2], and f{k) = otherwise. Here, e is a positive number, and 
\K\ is the width (the number of discrete wavenumbers) of the forcing band 
for the continuous (discrete) case. For this forcing, the spectral injection rate 
of ^a(^) is i^{k) = ^/\K\. Hence, the respective injection rates of ^'q, and ^'20 
are given by 
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jF{k)dk^e, (13) 

K 

j k'^F{k)dk^^ j k'^dk^s'^e, (14) 

where is the mean of k°' over the forcing band K. This forcing is used in the 
numerical simulations of Section 7, where e is a constant, yielding constant 
injection rates. 

Multiplying (6) by and subtracting the resulting equation from (7), noting 
that the forcing terms cancel, one obtains 



^(*2a - S"*a) =2l/^(s"*2M+, - *2M+2a)- (15) 

The difference on the right-hand side can be written as 



s"*2M+a - *2M+2a = -5'^(*2a - ^"^a) - / {s^^ - ^'^)(s" - k'^)^ ^{k) dk. 



Substituting this expression into (15) yields 



1 

'^2a - 5"*a) = -2i/^s2'^(^2a - s"^„) - j {s^>^ - k'"'){s" - k")<lf^{k) dk. 



dt 



For II > and a > 0, the integral in (16) is non-negative. The first order 
equation (16) is dissipative and forced by a non-positive term. It follows that 
^2a — s^^a becomes non-positive as i — > oo. Hence, the constraint 



*2a<s"*a (17) 

holds asymptotically. 

Remark 2. Actually, the constraint (17) holds for all t >T, where T is a finite 
time determined by the initial condition, except for some trivial cases, which 
are rather irrelevant in the context of turbulence (see [33]). This constraint 

is mathematically interesting in the bounded the inequality is in the 

opposite sense to Sobolev-type inequalities. In an unbounded domain where \E'q, 
may grow unbounded whether or not an inverse cascade exists, (17) becomes 
obvious. 



(16) 
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Remark 3. The dynamical constraint (17) holds if the single viscous operator 
is replaced by a combination A)^ for z/^, n > 0, where the summation 

is taken over all dissipation channels involved [33]. It also holds, with s replaced 
by k2, for a forcing satisfying k'^ F{k) dk < k^'Fik) dk < k^ F{k) dk. 
This double inequality automatically holds for any non-negative spectral in- 
jection, i.e. F{k) > for all k e K. 



4 Conservative transfer 

This section examines the nonlinear transfer of the quadratic invariants and 
their dynamical behaviours, particularly the possibility of a direct cascade 
of ^2a- A Fj0rtoft-type analysis [10] concerning the collective transfer of all 
nonlinear ly interacting triads is presented. The simple technique of Tran and 
Bowman [32] is employed in deriving a constraint on the spectral slope of the 
supposed direct-cascading range in terms of the growth rate of the inverse- 
cascading candidate and the dissipation wavenumber. These analyses finally 
lead to a self- consistent dynamical picture, which exhibits no direct cascade. 

4-1 Inverse vs. direct transfer 

Due to the conservation laws (5), the nonlinear transfer of and obeys 
certain restrictions. Suppose that an initial at wavenumber s, which corre- 
sponds to an initial = -^"^ai is nonlinearly redistributed in wavenumber 
space. Because of the conservation of '^2a, the redistribution of an amount ^+ 
beyond some given wavenumber k^ > s must satisfy < s°''^'^, or equiv- 

alently, ^^/^a < (•5/^+)"- Hence, no significant fractions of are allowed 
to be transferred toward wavenumbers k s. This condition corresponds to 
the well-known prohibition of a direct cascade of energy in the NS case. The 
constraint is more severe for a larger a in the sense that a larger a implies 
a smaller permissible ^'^ to be redistributed beyond a given /c+ > s. Simi- 
larly, due to the conservation of ^'q,, the redistribution of an amount 
to the wavenumbers lower than some given wavenumber k_ < s must satisfy 
s"\P^/A;" < \E'°, or equivalently, s"^~/(s"^°) < {k^/s)". Hence, no signif- 
icant fractions of s"^^ are allowed to be transferred toward wavenumbers 
k <^ s. On the other hand, the redistribution of any amount '^'^ < to- 
ward the vanishing wavenumbers is permitted. Such a transfer induces a loss 
of s"^'^ in ^2a- The conservation of then requires that the remainder 
\E'° — be redistributed toward the high wavenumbers, to make up for the 
lost amount s"^^. For a non-negligible remainder \1'^ — \1'^ and large a, the 
lost amount s"^'^ can be recovered with a redistribution of — "i/'^ toward 
wavenumbers slightly larger than s, resulting in no direct cascade. This means 
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that a direct cascade is possible only in the limit In other words, 

an inverse cascade can exist in the complete absence of a direct cascade, and 
the latter is not permitted until virtually all of ^'^ cascades to the vanishing 
wavenumbers. This observation is consistent with the results from numerical 
simulations of 2D NS turbulence that have found an inverse energy cascade, 
regardless of what happened to the enstrophy [4,5,12,32]. 



4-2 Constraint on spectral slope for k > s 



A more quantitative version of the preceding analysis is in order. The dual- 
cascade scenario is now examined in the framework of the forced-dissipative 
picture. A constraint on the spectral slope of the supposed direct-cascading 
range is derived in terms of the growth rate of and the dissipation wavenum- 
ber. 

The evolution of and is governed by 

2z^^*2^+a + / F{k)dk, (18) 
K 

2iy^^2^+2a + I k'^F{k)dk. (19) 

K 

Let us consider quasi-steady dynamics, where a steady spectrum has been es- 
tablished down to a wavenumber ko ^ s. The injections have become steady 
and satisfy F{k) dk = k"F{k) d/c, where s is a wavenumber in the forc- 
ing region. The evolution of has reached equilibrium, while continues 
to cascade toward wavenumbers /c < /cq at a steady growth rate d^a/dt — eq.^ 
It follows from (18) and (19) that 



Vtf z= — 

di " 



d ^ 



5%±2a^,a^^ (20) 

where e = Jj^ F{k) dk. For simplicity, it is assumed that the quasi-steady 
spectrum ^a(A;) can be approximated by [30] 



^ ^ iiko<k<s, , . 

^ It should be emphasized that althoTigh evidence exists of an inverse cascade for 
various values of a, one cannot be sure that an inverse cascade exists for all a > 0. 
Section 6 deals with the non-cascading case eg = 0, which is applicable to both 
bounded and unbounded systems in equilibrium. 
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where a, b, 7, (3 are constants, and ki, is the highest wavenumber in the range 
k^^, beyond which the spectrum is supposed to be steeper than k~^. The ratio 
on the left-hand side of (20) can be estimated as 

*2^+a - a /4 A;2M-7 dk + h /r k^^^-P dk 

_^„ Jws^^^^"'^dK + /,%,K^-^^-"d^ 

i/. '^''^-^ d/^ + /./^-^/^ d« ■ ^ ^ 

The inequahty results from dropping from the estimate of ^-2^+20 the spectral 
contribution beyond k^ and extending the range 6A;~^ to infinity in the estimate 
of The equality is obtained by making the respective substitutions 

K = k/s and k = s/k in the two integrals in each of the numerator and 
denominator plus the continuity relation as~'^ = bs~^. It follows that 



<— ■ (23) 



A direct cascade of requires that its dissipation occurs primarily in the 
small scales. Accordingly, the spectrum ^2^+2a(^") must be shallower than 
k~^. It is interesting to consider the regime where the onset of direct-cascading 
dynamics is possible, i.e., the regime for which ^2/i+2a(^) approximately scales 
as k~^ for k > s. This scaling means — /3 + 2iJ, + a ^ —1, or equivalently 
f3 — 2fi — 2^a — 1. The second integral in the denominator on the left-hand 
side of (23) is then of order unity (or smaller if a ^ 1). Now it is seen in 
Section 3 that a persistent inverse cascade (eo > 0) requires that ^'2//+a(^) 
be shallower than in the inverse-cascading range, i.e. 2/x — 7 > —1. This 
requirement implies that the first integral in the denominator on the left-hand 
side of (23) is of order unity. Hence, the denominator on the left-hand side of 
(23) is of order unity. On the other hand, the integral in the numerator can 
be significantly larger than unity. In particular, this integral approaches the 
value ln(A,v/s) in the limit /3 — 2// — a — 1 — > and can be large for a large ratio 
k,^/s. This result implies that the left-hand side of (23) can be significantly 
larger than unity in the non-dircct-cascading regime. The constraint (23) then 
requires that eo be sufficiently close to e. Hence, an inverse cascade with various 
strength (various values of the ratio eo / e) can exist in the absence of a direct 
cascade, as was previously argued in this section without an explicit dissipation 
mechanism. Finally, in the limit s/A;,^ — > 0, as is commonly assumed for high- 
Reynolds number turbulence, a direct cascade, regardless of the spectral slope 
— /3, requires the limit eo — > e. 

Remcirk 4. For NS turbulence, the critical slope of the energy spectrum. 
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for /c > s, at the onset of a possible direct cascade is —5. The KLB direct- 
enstrophy-cascading range scales as A';^'\ which is mnch shallower than the 
critical scaling k~^. (The spectral discrepancy between these two scaling laws 
is huge.) /? has been proposed to have some favourable values between 3 and 5, 
according to the theories of Moffatt [20] (/3 = 11/3), Saffman [24] (/3 = 4), and 
Sulem and Prisch [29] {P < 11/3). For SQG turbulence, the critical slope is 
—2. The Kolmogorov-Kraichnan direct-energy-cascading range scales as k~^^^, 
which is slightly shallower than the critical scaling k~'^. 

The question remains as to whether a direct cascade is realizable in unbounded 
fluids. The present analysis narrows down the possibility of a direct cascade 
to the regime in which the growth rate of \1/q, approaches its own injection 
rate. If a direct cascade could be realizable and if a spectrum much shallower 
than the critical spectrum could be achieved (e.g., a k'^ scaling in NS turbu- 
lence) , the critical scaling would intuitively be realizable at relatively moderate 
Reynolds numbers, and the direct-cascading range would become shallower as 
the Reynolds number increases. Numerical simulations may be able to probe 
into this critical regime and help assess the existence or non-existence of a 
direct cascade. Given the low resolutions available in current computers, a 
plausible possibility is to explore how eo, /3, and ki, adjust with respect to the 
Reynolds number in the non-direct-cascading regime [32]. This information 
may then be extrapolated to the critical limit. A difficulty with this scheme 
is that one might not be able to get sufficiently close to the critical regime to 
make the extrapolation meaningful. ^ 



4-3 A non- direct- cascading dynamical picture 

It is interesting to consider a self-consistent dynamical picture, which is in 
accord with the preceding analysis and exhibits no direct cascade, even in 
the limit of infinite Reynolds number. This discussion is concerned with the 
plausibility of a direct cascade rather than with an attempt to find an alter- 
native to the dual-cascade theory. In the limit i/^ — > 0, the equilibrium of 
implies that '^2fi+2a diverges as In the absence of a direct cascade, 
this condition can be accounted for by both the logarithmic divergence of 
ln(Av/s). as A;^ oo. and an increase of b. The latter factor is relevant since 
it is not expected that k,^ grows exponentially as z/^ 0. In any case, one has 
s"e/2i/^ = 'if2n+2a > bln{k^/s). Meanwhile, is given by 



^ For 2D NS turbulence, one would need several wavenumber decades to simulate 

an inverse energy cascade that carries about 90% of the injected energy to the large 
scales [32]. [This requirement can be readily deduced from (25).] This resolution is 
already a challenge for current computers. In fact, no numerical simulations per- 
formed thus far have achieved a resolution better than four wavenumber decades. 
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= hs-^ —— + - , (24) 

\2/i + l-7 a) 

for 2/i + 1 — 7 > 0. The dissipation of can be bounded from above by 




This result imphes that 2i/^^'2^+q — > as — > oo (i/^ — > 0). Hence, in 
the hmit — > 0, no direct cascade would be needed in order for an inverse 
cascade to proceed toward A: = 0, carrying with it virtually all of the input of 
Thus, this hypothetical picture depicts non-direct-cascading dynamics for 
all Reynolds numbers. As becomes ever smaller, the range bk~^ gets ever 
wider, (3 gets ever closer to the critical value (3 = 2ii + a + 1, and h increases 
somewhat less rapidly than Eventually, (for all real 9) diverges as 
— > 0. This feature is completely absent from the dual-cascade picture, ^ 
where \E'2^+a is supposed to remain finite in the limit 0. Note that 

the present discussion would apply to both bounded and unbounded systems. 
The only essential difference between these systems would be that for the 
bounded system, the inverse cascade is halted at the largest available scale, 
upon which significant departures in the dynamics of the two systems take 
place. The approach to equilibrium after the arrest of the inverse cascade 
in 2D NS turbulence was numerically studied by Borue [5]. Also, a physical 
interpretation of the observed k^^ spectrum at the large scales was given. Tran 
and Bowman [30,32] observed and explained a similar result. 



The existence of a weak inverse energy cascade in the complete absence of 
a direct enstrophy cascade for moderate-Reynolds number NS turbulence is 
strongly supported by the theoretical arguments and numerical results of Tran 
and Bowman [30,32] and by the numerous numerical investigations targeted 
at the inverse cascade in the literature. ^ Section 7 provides further numerical 
evidence for the realizability of non-direct-cascading dynamics. The results 
reported there also show the formation of an enstrophy range shallower than 
k"^ during the transient dynamics, thereby suggesting that the turbulence 



^ In the dual-cascade picture, b is essentially insensitive to in the limit 0. 
The divergence of ^2n+2a, as 0, is mainly due to the increase of ky. (and 
possibly other quadratic quantities) remains finite in the limit ^ 0. Therefore, 
the dual-cascade picture poses a "sudden" blow-up at = 0. This behaviour is 
somewhat harder to perceive than the "smooth" blow-up in the non-direct-cascading 
picture. 

^ Those investigations inadvertently do not recognize the absence of a direct cascade 
for obvious reason: the forcing region is usually chosen as close as possible to the 
largest available wavenumber to obtain a wide inverse-cascading range, thereby 
allowing for no simultaneous observations of the enstrophy dynamics. 
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would become non-direct-cascading by evolving through a direct-cascading 
stage. 



5 Non-conservative transfer 



The creation and depletion of quadratic quantities other than the two in- 
variants by nonlinear transfer are now investigated. For an arbitrary spec- 
trum, subject to the finiteness of the quadratic quantities under considera- 
tion, one has the following interpolation-type inequality, which can be shown 
by Holder's inequality: 



<\+e,<<\+e,<\-'\ (26) 

for ^1 > ^2 > and 9^ is an arbitrary real number. This inequality is geometry- 
independent, so that it applies to both bounded and unbounded cases. The 
equality occurs if and only if \t'(A;) 7^ is monoscale, i.e., if ^(/c) is nonzero only 
at a single wavenumber (a delta function for the case of continuous spectrum). 
This trivial case is excluded so that only the strict version of (26) is applied in 
what is to follow. By applying (26) with (^1, 82, Os) — {a,9 — a, a) one obtains 



*e<*2a''"'*a"'^" for ee{a,2a). (27) 

Similarly, the two substitutions (^1,^2,^3) = (^ ~ a, a, a) and {Oi, 02,03) — 
{2a -9,a-9,9) into (26) yield 



^g>^Z''-'^l-'/'^ for e^[a,2a]- (28) 

The quantities on the right-hand sides of (27) and (28) are conserved by non- 
linear transfer (because ^'^ and ^^20 are invariant) and equal to ^'g, for an 
initial monoscale ^P^. Therefore, the subsequent nonlinear transfer necessarily 
results in the depletion (creation) of ^0 for 9 e (ct, 2a) {9 ^ [a, 2a]). It should 
be noted that continuous increase (decrease) of the non-conserved quadratic 
quantities from an initial monoscale spectrum is possible but cannot be as- 
sured, due to the reversible nature of the inviscid problem. 

The above picture, when applied to forced-dissipative systems in equilib- 
rium or in quasi-steady dynamics, implies that the nonlinear transfer de- 
stroys (enhances) the injection of \E'e from a monoscale forcing for 9 G {a, 2a) 
{9 ^ [a, 2a]). In other words, the integrated transfer k^^°'Ta{k) dk in the 
evolution equation 
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d_ 
dt 



oo 

^9 = j k'-^Ta{k) dk - 211^2^.^9 + J k'-''F{k) dk (29) 

K 



is negative (positive) for 6 G {a.. 2a) {9 ^ [a, 2a]). To see this result explicitly, 
note that in equilibrium or quasi-steady dynamics, Ta{k) > for 7^ s because 
F(k) — ioi k s, and the dissipation is non-negative [of. (4)]. Similar to 
(26), the non-negative (and non-monoscale) transfer function Ta{k), for k ^ s, 
satisfies 

I k'^'+'^'Ta{k)dk\ < j k''+'''Ta{k)dk\ I k''Ta{k)dk\ (30) 

for ^1 > > and ^3 is an arbitrary real number. Upon substituting 
(^^1, ^^2, ^'3) = {a,9 — a,0) into (30) and using the conservation laws (5), one 
finds 



J k^-'^Taik) dk < 

k^s 



0/a-l / \ 2-0/a 



k^'Taik) dk 



Ta{k)dk 



^is^\T^is)\f'^-' \T^is)\'-'/^ = s'-^\T^is)\ (31) 
for 6* e (a, 2a), or equivalently, 

00 

J k''-''Ta{k)dk<0. (32) 


This result means that for 9 G {a, 2a), ^'e* is destroyed by nonlinear transfer. 
Similarly, with the substitutions (^1,^2,^3) — {9 — a,a,0) and (^1,^2,^3) = 
{2a — 9, a — 9,9 — a), one obtains 

00 

J k''-''Ta{k)dk>0, (33) 


for 9 ^ [a, 2a]. This result means that for 9 ^ [a, 2a], is created by non- 
linear transfer. 

The creation and depletion of for different values of a are interesting; here, 

only some typical cases are considered. For a < 2, the enstrophy ^4 (and 
vorticity) can be created. This enstrophy creation arguably occurs in SQG 
turbulence (a = 1), as its inviscid unforced dynamics is known to resemble 
that of the three-dimensional Euler equation in many aspects, particularly 
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in the possibility of spontaneous development of singularities [7,9]. Note that 
the enstrophy creation need not be a consequence of a direct cascade and 
that the created enstrophy is subject to viscous dissipation in dissipative sys- 
tems. For 2 < a < 4 such as in RSF turbulence {a — 3), the energy is created 
while the enstrophy is destroyed. The creation of energy occurs toward the low 
wavenumbers and evades viscous dissipation. If an inverse cascade of ^'q, pro- 
ceeds toward A; = at a steady growth rate, the energy growth rate increases 
without bound. When a > 4, both the energy and enstrophy are created, and 
this creation occurs toward the low wavenumbers. Both the energy and enstro- 
phy growth rates increase without bound in the presence of a steady growth 
rate of "^a- When a < 1, both the energy and enstrophy are created, and 
this creation occurs toward the high wavenumbers. The NS system happens 
to have the critical value o; = 2, for which both the energy and enstrophy are 
conserved by advective nonlinearities. 



6 Constraints on the slopes of equilibrium spectrum 



For fluids in equilibrium, eg = and (20) reduces to the balance equation 
s°"^2ti+a = ^2^+2o- Assuming the spectral scaling (21) and following the cal- 
culations leading to (23), one obtains the estimate of '^f2^l+2a — s"'if2ij,+a- 



J {k" - s")A;2''^„(A;) dk>aj{k''- s'')k^>'-^ dk + b J {k"" - s'')k^^-^ dk (34) 

s 

1 



ko ko 

1 



= as 



ko/s 

1 

s/fc,y 



J (l-0«'''~''d«+ / {1 - k'')k^-^^-''-^ dK 

\ ko/s s/kv , 



where the inequality results from dropping the spectral contribution beyond k^,, 
the second line is obtained by the respective changes of variables K — k/s and 
K — s/k in the two integrals on the right-hand-side of the first line, and the 
third line is obtained by using the continuity relation as~^ = bs~^ . Here, 
is the lowest wavenumber in the bounded case or a cutoff wavenumber in the 
unbounded case. Since the left-hand side of (34) vanishes, it can be deduced 
that 
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1 1 



{1- K'^)K^''-^dK> / (l-O 



■P-^f'-'^-^ dK. 



(35) 




For bounded high-Reynolds number turbulence driven at a relatively low 
wavenumber s, it is reasonable to assume that kg/s > s/ku. This condition 
requires 2/x — 7 < (3 — 2/1 — a — 2, as the integrals in (35) decrease if the 
corresponding powers of k, (2/x — 7 and P — 2/1 — a — 2) increase. It follows 
that 7 + /3>4// + q; + 2. 

On the other hand, the convergence of the right-hand integral in (35), as 
s/ki, — >• 0, requires P > 2/1 + a + 1. (Equivalently, ^'2/i+2a(^) is steeper than 
for k > s.) That is, no direct cascade exists for equilibrium dynamics. The 
spectral contribution beyond ki, which is dropped from (34) is then negligible 
and (35) is essentially an equality. It follows that 'j + (3 ~ 4^ + a + 2 if 
ko/s ^ s/ki,, and 7 + /9 can fall slightly below 4fi + a + 2 if k^/s < s/k^. 
For both these cases, 7 < 2/i -|- 1, so that an inverse cascade is also ruled out 
because the dissipation of cannot occur at the largest scales. Although an 
inverse cascade cannot be ruled out in the case k^/ s > s/k^,, such a cascade, 
if realizable, would only be marginal and characteristically different from its 
unbounded counterpart in the classical picture. 

It is more informative to write the constraint 7 -|- /3 > 4// -|- a -|- 2, for kg/s > 
s/ky, in terms of the slopes of the kinetic energy spectrum ^'2(^)- Let —7 and 

—{3 denote, respectively, the slopes of ^2(^) in the ranges of wavenumbers 
lower and higher than the forcing wavenumber; 7 and j3 are then related to 7 
and /3by7 = 7-|-a — 2 and 13 = (3 + a — 2. \i follows that 



This inequahty reduces to 7 -|- /5 > 8 for NS turbulence, 7-|-/9>3forSQG 
turbulence, and 7 -|- /9 > 4// -|- 7 for RSF turbulence. 

Constantin [8] has studied the spectral scaling of SQG turbulence and derived 
the pointwise estimate '^2{k) < ck~^ in the wavenumber range k < s. Smith 
et al. [28] numerically explored this problem and found energy spectra con- 
siderably shallower than A;^^ in this range. The "energy" spectrum plotted by 
Smith et al. [28] is ^i(/i;), and this spectrum appears to have a slope ^ —1.5. 
Hence, ^2(A;) approximately scales^ as k^^-^ . Dimensional analyses predict the 
scaling ^2{k;) oc k~^ ioi k < s [23,28]. The present estimate agrees with Smith 

^ Note that Smith et al. [28] used Unear drag in their simulations. Hence, the 
observed k^^'^ scaling may in principle differ from that of SQG turbulence with the 
dissipation operator (— A)^/^. Nevertheless, their results provide strong evidence 
that the energy spectrum in the wavenumber range A; < s is shallower than k"^. 



7 + /? > 4/x + 3q; - 2. 



(36) 
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et al.'s numerical results and with the prediction based on dimensional analy- 
ses better than with Constantin's estimate. On physical grounds, a spectrum 
shallower than is more plausible, for two reasons. First, the prohibition of 
a significant inverse energy transfer in SQG turbulence (in marked contrast to 
the inverse energy cascade in NS turbulence) should render the energy spec- 
trum in the wavenumber range k < s shallower than because any spectra 
not shallower than k~^ would require a significant inverse energy transfer. 
Second, the hypoviscous dissipation of SQG dynamics, (— A)^/^, does not ac- 
commodate any energy spectra not shallower than k~^ in the low-wavenumber 
range (see Section 3). 

In a recent laboratory experiment on a rapidly rotating fluid, Baroud et al. [1] 
observed a k^"^ spectrum extending for almost two wavenumber decades lower 
than the forcing wavenumber. This result is consistent with the constraint 
7+/3 > 8 obtained for the NS case, but cannot be explained within the context 
of SQG turbulence, for the reasons discussed above. A possible explanation 
for Baroud et al.'s result is that a rotating fluid, believed to be governed by 
the SQG equation, ® may still behave (to some extent) as a 2D NS fluid. 
Therefore, the dynamics of such a fluid might exhibit mixed characteristics of 
both SQG and NS turbulence. While a systematic analysis is required for a 
proper understanding of Baroud et al.'s experimental result, this hypothesis 
offers a possible explanation for the observed spectrum. 



7 Numerical results 

This section presents results from numerical simulations that confirm the re- 
alization of an inverse energy cascade in 2D NS turbulence, where energy is 
transferred to the large scales via the Kolmogorov-Kraichnan spectrum k~^^^, 
in the absence of a direct cascade. These results provide numerical evidence 
(also see [32]) for the prediction and motivation for the discussion in the pre- 
ceding sections. For the two sufficiently resolved Reynolds numbers, the energy 
growth rate cq, calculated as the inverse cascade reaches k ^ 3,is larger for the 
higher Reynolds number. However, this preliminary attempt is in no position 
to address what happens in the high-Reynolds number regime. 

Eq. (1) with a — 2 and /i — 1 is simulated in a doubly periodic square of side 

^ The derivation of the SQG equation involves the temperature field [22,34]. More 
accurately, in order to arrive at the SQG equation, the vertical gradient of the 
temperature at the surface of a rotating stratified fluid is considered. Such a gradient 
is absent in Baroud et al.'s experiment since the fluid temperature is presumably 
constant throughout the rotating tank. Hence, the results in [1] cannot be attributed 
with certainty to SQG turbulence. 
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27r. The forcing f{k) is given by (12), where e = 1 and K = [49.5,50.5]. The 
enstrophy injection rate is s^e ~ 2500, where ^ 2500 is the mean of k"^ over 
K. The attractive aspect of (12) is that it is steady, so that the uncertainty 
in the energy transfer rate eo, which is calculated according to eo = e — 21/1^^4, 
is due solely to that of the enstrophy ^4. It is necessary to have a sufficiently 
wide inverse-cascading range, so that the enstrophy approaches its equihb- 
rium value before finite-size effects interfere with the inverse cascade. It is 
also desirable to have a wide wavenumber range A; > s, so that high Reynolds 
numbers can be achieved. For a given resolution, the choice of is a com- 
promise between these two considerations. With the present choice of K, the 
enstrophy calculated for a k~^^^ spectrum in the wavenumber interval [3, 50] 
is approximately 2% less than that calculated for the same spectrum in the 
interval [0,50]. Hence, the energy transfer rate, which is calculated when the 
energy peak reaches wavenumber /c ?a 3, may be greater than its projected 
value by a few percents (assuming no significant finite-size interference to the 
inverse cascade up to that point in time). Dealiased pseudospectral method is 
used with 1024^ and 2048^ total modes. Several different Reynolds numbers, 
corresponding to several different values of Ui, were simulated. The lower reso- 
lution was used for i/i = 7 x 10""^; 5 x 10""^, and the higher resolution was used 
for i/i = 3.5 X 10-^; 2.5 x 10"^; 1.25 x 10"^; 6.25 x 10"^ For all simulations, the 
initial energy spectrum had a negligibly small peak at the forcing wavenumber 
(*2(50) ^ 10-1°). 

For ui = 7 X lO^'i, the energy peak was observed to reach /c 3 at t ~ 10.2, 
upon which both the energy spectrum (see Fig. 1) and the total enstrophy 
were calculated (both time-avcragcd between t = 10.2 and t — 10.3). For a 
sense of this time scale, the dissipation time at the forcing wavenumber s ^ 50 
is ts = {2uis'^)^^ ~ 0.3. The calculated value of enstrophy \E'4 = 588 amounts 
to an energy dissipation rate 2t/i^'4 = 0.823, accounting for 82% of the energy 
injection rate e — 1. The energy growth rate eo is then 18% of the energy 
injection rate, due to the inverse energy cascade. It is evident from Fig. 1 that 
an inverse cascade is realizable in the complete absence of a direct cascade: 
the enstrophy range scales as k~^''^, so that virtually all of the enstrophy is 
dissipated in the vicinity of the forcing region. For z/i = 5 x 10"^, the energy 
peak was observed to reach A; ?a 3 at t 8.8, upon which both the energy 
spectrum (see Fig. 2) and the total enstrophy were calculated (both time- 
averaged between t = 8.8 and t = 8.9). The calculated value of enstrophy ^4 = 
761 amounts to an energy dissipation rate 2z/i^4 = 0.761, accounting for 76% 
of the energy injection rate e = 1. The energy growth rate eo is then 24% of the 
energy injection rate, due to the inverse energy cascade. The enstrophy range 
scales as k'~^^^, which is slightly shallower than its counterpart in the previous 
case. These results imply that the inverse cascade and the enstrophy-range 
spectrum get stronger and shallower, respectively, as the viscosity coefficient 
z/i decreases (the Reynolds number increases). Finally, for both cases, the 
inverse cascade carries only a small fraction of the energy input to the largest 
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scale and yet a k~^/^ spectrum manifests itself nonetheless. This suggests that 
a k~^/^ inverse-cascading range does not require the transfer of virtually all 
energy input to the large scales. 




10" 10' 102 

k 

Fig. 1. The energy spectrum *2(fc) vs. k averaged between t = 10.2 and t = 10.3. 

The cnstrophy averaged in the same period is ^'4 = 588. The energy injection rate 
and viscosity coefficient used are e = 1 and z^i = 7 x 10^^, respectively. The inverse 
energy transfer rate is eo = e — 2z/i^4 = 0.823, which accounts for 82% of the energy 
injection rate. The enstrophy range scales as k"^-"^. Hence, an inverse energy cascade 
carrying 18% of the energy injection rate is realized in the complete absence of a 
direct cascade. The theoretical predictions are discussed in Section 4. 

The steep enstrophy-range spectra in Figs. 1 and 2 seem to suggest that 
the present resolutions, particularly the higher one with 2048^ total modes, 
may be sufficient to resolve higher Reynolds numbers. It turns out, how- 
ever, that further decrease of ui leads to spurious dynamics. For example, 
for Ui = 3.5 X 10^^; 2.5 x 10^^, the kr^^^ inverse-cascading range is distorted 
after the energy peak has reached k ^ 8 {at t ^ 3.5): the energy range be- 
comes shallower, due to a spectral rise at the forcing region. It is not known 
whether this phenomenon is due to finite-size or finite-resolution effects or 
both. (Such a rise is not observable in the previous simulations, where an 
inverse-cascading range is established down to A; ~ 3 without any signifi- 
cant departure from the expected —5/3 slope.) A reason to suggest that this 
phenomenon is mainly a manifestation of finite-resolution effects is that in 
these cases, the palinstrophy quickly grows to exceed its equilibrium value 
^6 = s^e/2i/i with most of the (resolved) enstrophy range shallower than k"^ 
(see Fig. 4, for example), indicating insufficient resolution for the transient 
dynamics: a significant amount of palinstrophy would be carried by the unre- 
solved (truncated) wavenumbers. The subsequent evolution of the enstrophy 
range, confined within the insufficient-resolution range, was observed to cause 
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Fig. 2. The energy spectrum ^'2(fc) vs. k averaged between t = 8.8 and t = 8.9. 
The enstrophy averaged in the same period is *4 = 761. The energy injection rate 
and viscosity coefficient used are e = 1 and vi = 5 x 10"^, respectively. The inverse 
energy transfer rate is eo = e — 2ui^!/^ = 0.761, which accounts for 76% of the energy 
injection rate. The enstrophy range scales as k~^'^. Hence, an inverse energy cascade 
carrying 24% of the energy injection rate is realized in the complete absence of a 
direct cascade. The theoretical predictions are discussed in Section 4. 

a rise of the spectrum around s. A growth of enstrophy is induced, thereby 
interfering with the inverse energy cascade. This phenomenon may be viewed 
as a redistribution of a significant palinstrophy amount back to the forcing re- 
gion. This dynamical behaviour is illustrated by Fig. 3, from which the energy 
range is seen to deviate considerably from the Kolmogorov-Kraichnan k~^/^ 
spectrum. 

Higher Reynolds numbers seem to produce shallower enstrophy range during 
the transient dynamics. For Vi = 1.25 x lO^'^ (z/i = 6.25 x 10^^), Fig. 4 
(Fig. 5) shows the energy spectrum averaged between t = 1.1 (t = 0.90) and 
t — 1.2 {t = 0.95), during which time about 0.7ts (0.3ts) has elapsed. The 
transient enstrophy range scales as A;"^'^ {k~^-^) and extends to k^, ^ 500 = 10s 
(ki, > 10s). ^ The enstrophy dynamics have become fully dissipative since the 
palinstrophy has already reached \E'6(1-1) ~ 1-05 x lO'^ (^6(0-9) ^ 2.15 x lO'^), 
slightly exceeding its equilibrium value s^e/2vi = 1 x 10^ (s^e/2z/i = 2 x 10^). 
Clearly, the present resolution is insufficient, and the subsequent evolution is 



^ Since the formation of an enstrophy range steeper than (but possibly arbitrar- 
ily close to) k~^ requires virtually no enstrophy to be transferred to its higher- 
wavenumbcr end, the enstrophy range may acquire a slope between —5 and —3 
during the transient stage. The width and slope of this range should depend on the 
Reynolds number and evolve with time. 
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Fig. 3. The energy spectrum ^'2(fc) vs. k averaged between t = 7.74 and t = 7.84. 
The enstrophy injection rate and viscosity coefficient used are s^e 2500 and 
ui = 2.5 X 10~^, respectively. The energy range is significantly shallower than 
k^^^^, due to a rise of the spectrum around s. This rise occurs (and persists) as 
the palinstrophy evolves within the insufficient-resolution range, after exceeding its 
equilibrium value 5 x 10^ (achieving the value 5.26 x 10^ at t ^ 1.2) with an 
enstrophy-range slope —5. 



supposed to suffer from what have been termed "finite-resolution effects" as is 
seen in the case ui — 2.5 x 10"^. Apparently, the non-direct-cascading picture 
would not apply to this stage, nor would it apply to the subsequent dynamics 
if a persistent direct cascade could be established, withstanding any significant 
adjustments of the cnstrophy-range spectrum. Rather, the realization of this 
picture requires that the enstrophy range relax toward a k^^ spectrum. This 
relaxation would essentially involve a redistribution of palinstrophy (which 
has already achieved its equilibrium value \l/6 ^ s^e/2z/i) in the wavenumber 
range k > s. In some sense, the realizability of the non-direct-cascading picture 
requires that any direct-cascading states, regardless of the enstrophy-range 
slope, be "unstable" and that the palinstrophy be "free" to redistribute itself, 
until a spectrum can be achieved. This dynamical description is supposed 
to be what subsequently occurs to the transient states depicted by Figs. 4 
and 5 (provided sufficient resolutions). In any case, the realization of transient 
direct-cascading dynamics, for ui = 2.5 x 10^^; 6.25 x 10^"'', suggests that for 
high Reynolds numbers, a transient direct-cascading stage would precede the 
non-direct-cascading picture if such a picture could be realizable. 
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Fig. 4. The energy spectrum ^'2(fc) vs. k averaged between t = 1.1 and t = 1.2. The 
palinstrophy averaged in the same period is 1.05 x 10^. The respective enstrophy 
injection rate and viscosity coefficient used are s^e ^ 2500 and ui = 1.25 x 10~^, 
which amount to the pahnstrophy equihbrium vahic 1 x 10^. One observes that the 
palinstrophy exceeds its equilibrium value well before a spectrum is formed 

in the energy range. The enstrophy range is shallower than k~^, clearly indicating 
insufficient resolution for the subsequent evolution. 

8 Concluding remarks 

The nonlinear transfer and asymptotic dynamical behaviours of a turbu- 
lence have been investigated in this work. The results obtained include sev- 
eral estimates on the dynamical quantities, such as the dissipation of the 
inverse-cascading candidate, and various constraints on the spectral distri- 
bution of energy. The conservative transfer of quadratic invariants and the 
non-conservative transfer of other quadratic quantities are investigated. The 
possibility of a direct cascade in unbounded fluids is examined by exploring 
the constraint imposed by the conservative transfer of the advective term. 
A condition for the possible onset of direct cascade is derived, from which 
the dual-cascade hypothesis is reformulated. The realization of an inverse en- 
ergy cascade in the complete absence of a direct enstrophy cascade in 2D NS 
turbulence is confirmed by numerical simulations. 

For systems driven by a forcing spectrally localized at wavenumber s (in the 
sense discussed in Section 3) and dissipated by a viscous operator A)'*, 

where z/^ > and > 0, the invariants and satisfy < 
The nonlinear transfer destroys (creates) the quadratic quantity for 6 G 
{a, 2a) {6 ^ [q;,2q;]). This phenomenon allows for vorticity to be created 
for systems with a < 2, which include the SQG equation. The enstrophy 
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Fig. 5. The energy spectrum ^'2(fc) vs. k averaged between t = 0.90 and t = 0.95. 
The pahnstrophy averaged in the same period is 2. 15 x 10^. The respective enstro- 
phy injection rate and viscosity coefficient used are s^e ^ 2500 and ui = 6.25 x 10~^, 
which amount to the pahnstrophy equihbrium vahic 2 x lO'''. One observes that the 
pahnstrophy exceeds its equihbrium value weU before a spectrum is formed 

in the energy range. The enstrophy range is shallower than k~^, clearly indicating 
insufficient resolution for the subsequent evolution. 

creation occurs in the high wavenumbers, and the created enstrophy is subject 
to viscous dissipation. For systems with a > 2, such as the RSF equation, 
energy is created. The creation of energy in RSF turbulence occurs in the 
low-wavenumber region. In the presence of a persistent inverse cascade, this 
energy creation would occur at an ever-growing rate. 

Fj0rtoft-type arguments, based on the two conservation laws of 2D NS turbu- 
lence, are reformulated for a turbulence, in an attempt to gain more insights 
into the possibility of a direct cascade. It is shown that a weak inverse cas- 
cade, one that does not carry virtually all of the injection of ^q, to ever-larger 
scales, cannot co-exist with a direct cascade of "if 2a- In order for a direct 
cascade to be permitted (not to say that it actually occurs), the inverse cas- 
cade is required to proceed toward the vanishing wavenumbers, carrying with 
it virtually all of the injection of This strong inverse cascade is widely 
believed to be what happens in the limit of infinite Reynolds number, ac- 
cording to the classical dual-cascade hypothesis. However, it is argued in this 
present paper that a dynamical picture that features no direct cascade for all 
Reynolds numbers is self-consistent and possible. This result considerably nar- 
rows down the possibility of a direct cascade and also explains the robustness 
of the inverse energy cascade observed in numerical simulations of 2D NS tur- 
bulence, regardless of what happens to the enstrophy [4,5,12,32]. Furthermore, 
this hypothetical picture serves as a common dynamical background for both 
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bounded and unbounded fluids. In both cases, no direct cascade would occur 
for all Reynolds numbers, and differences would arise only when the inverse 
cascade gets reflected by the largest available scale in the bounded case. 

For a turbulence confined to a periodic domain, with a viscous dissipation 
operator (— A)'^ (/j, > 0) the only dissipation mechanism, no direct cascade 
is permitted. Moreover, an inverse cascade is not possible if the range of 
wavenumbers lower than the forcing wavenumber is sufficiently wide. More 
precisely, if the forcing wavenumber s is no smaller than the geometric mean 
^/koh> of the integral wavenumber ko (corresponding to the system size) and 
the dissipation wavenumber ki, (the wavenumber beyond which the spectrum 
becomes steeper), then an inverse cascade is not possible. For s < ^/koky, an 
inverse cascade cannot be ruled out; however, such a cascade would only be 
marginal and characteristically different from the inverse cascade in the classi- 
cal picture for unbounded systems. These results are the generalized versions 
of those obtained in [30,33] for the NS system. They also apply to unbounded 
fluids in equilibrium (if such an equilibrium could be achieved), where k^ is a 
wavenumber cutoff. 

Numerical analyses on 2D NS turbulence have been performed to conflrm that 
a weak inverse energy cascade is realizable in the complete absence of a direct 
enstrophy cascade. In particular, it is shown that an inverse energy cascade, 
carrying approximately 24% (18%) of the energy injection to the large scales, 
is realized and accompanied by a k~^'^ {k~^^'^) enstrophy range. The results 
also suggest that the inverse- cascading strength (enstrophy-range spectrum) 
becomes stronger (shallower) as the Reynolds number increases and that for 
high Reynolds numbers, the non-direct-cascading picture (if it could be realiz- 
able) would be preceded by a transient direct-cascading stage. Unfortunately, 
due to resolution limitation, the numerics fall short on addressing quantita- 
tively how the non-direct-cascading picture evolves with progressively higher 
Reynolds numbers. This open problem undoubtedly holds the key to a better 
understanding of the many aspects of flnite- Reynolds number turbulence, pos- 
sibly including the dual-cascade hypothesis in the limit of inflnite Reynolds 
number. In the most favourable scenario, in which the KLB theory should 
turn out to be correct, a knowledge of this problem would be an asset for for- 
mulating a more complete theory of 2D NS turbulence. In particular, the new 
theory would have to address the dynamics corresponding to the (huge) spec- 
tral "gap" between the critical —5 power-law scaling and the —3 power-law 
scaling of KLB theory (with or without a logarithmic correction) . Intuitively, 
this "gap" would have to be fllled with progressively higher Reynolds num- 
bers, starting from the critical Reynolds number corresponding to the onset 
of a direct enstrophy cascade. 
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